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ABSTRACT 


Conventionally, the control of liquid slosh system is done based on model- 
based techniques that challenging to implement practically because of the 
chaotic motion of fluid in the container. The aim of this article is to develop 
the tuning technique for model-free PID with derivative filter (PIDF) 
parameters for liquid slosh suppression system based on particle swarm 


optimization (PSO). PSO algorithm is responsible to find the optimal values 
for PIDF parameters based on fitness functions which are Sum Squared Error 
Keywords: (SSE) and Sum Absolute Error (SAE) of the cart position and liquid slosh 
angle response. The modelling of liquid slosh in lateral movement is 
considered to justify the design of control scheme. The PSO tuning method is 
compared by heuristic tuning method in order to show the effectiveness of 
PID controller the proposed tuning approach. The performance evaluations of the proposed 
PSO tuning method are based on the ability of the tank to follow the input in 
Sum squared error horizontal motion and liquid slosh level reduction in time domain. Based on 
the simulation results, the suggested tuning method is capable to reduce the 
liquid slosh level in the same time produces fast input tracking of the tank 
without precisely model the chaotic motion of the fluid. 
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1. INTRODUCTION 

The movement of the container system which contains liquid inside it usually faced a problem due 
to motion of the free liquid surface as known as slosh. The main problem of liquid slosh is it produces an 
additional moments and forces which can disturb the performance of the system. Hence, liquid slosh 
suppression is essential to solve the problems in numerous areas. For instance, the dynamic behavior of the 
ship carried a partly liquid filled tank onboard, is always disturbed because of additional forces and moments 
produced by a liquid slosh [1]. Another example, the pouring work in metal industries usually done by 
human operator. As we know, human tend to do an error because of many factors such as tiredness and 
carelessness. Consequently, human can cause sloshing and the molten metal will spill out from the ladle [2]. 
Moreover, the dynamics of the moving liquid cargo also can be troubled by the liquid slosh forces and 
moments acted in vertical and horizontal directions [3]. Thus, it is vital to suppress this liquid slosh problem 
during the movement of the container. 

The goal of this study is to design the controller for the liquid slosh suppression HIL system so the 
tank can follow the prescribed location precisely in the same time, minimize the slosh angle. In recent years, 
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several studies have focused on the liquid slosh control. Due to constraints of liquid slosh measurement, 
several passive approaches are executed to control the liquid slosh. For instance, absorbers and baffles are 
used as a passive elements to reduce the slosh energy in the container [4-5]. The drawbacks of these 
approaches to the system are heavy, bulky, require much time to construct and complex. Several attempts 
have been made to suppress the slosh of the liquid in the container by using feed-forward controller. 
For instance, minimum time feedforward control [6], filtering techniques [7], hybrid command smoothing 
and input shaping [8] and input shaping technique [9, 10]. These methods are used by generating the 
prescribed motion without using any feedback sensors, which reduced the residual slosh of the liquid 
generated in the tank. Disappointingly, feed-forward controllers are very sensitive to external perturbation 
occurred in the system. As an alternative, closed-loop control or feed-back control, which is popular with its 
robustness to external perturbation, has been applied to reduce the slosh of the liquid. For example, Variable 
Gain Super-twisting Algorithm (VGSTA) [11], sliding mode control [12], PID control [13-15], H, 


control [16], active force control (AFC) [17] and single input fuzzy logic controller [18]. 

As we know, practically, the traditional control techniques which depends on precise modelling of 
the system as known as model-based control, are hard to implement in the liquid container system because of 
difficulties to model the dynamic motion of the liquid and turbulent behavior of liquid slosh in the container. 
So, the more attractive approach which the controller design is done without knowing the precise model of 
the system as known as model-free control, has a good potential to implement in liquid container system for 
reducing the slosh. One of the model-free controller is Proportional-integral-derivative (PID) controller. 
PID is well known with its reliability and robustness and it has been widely used on the numbers of 
applications such as in [19-21]. Despite PID controller is able to solve those system control problems, 
the tuning of PID parameters is very tedious task. So, it is important to decrease the tuning time for PID 
controller to take full benefits of it. Therefore, for the last few decades, the optimal tuning of PID parameters 
are done by implementing the optimization techniques that proven to reduce the time of PID tuning. 

The main objective in this paper is to find the optimal value of PIDF controller parameters for the 
nonlinear liquid slosh model. PSO is chosen to tune 8 parameters; Kp, Ki, Ka and N to the optimal values for 
two PIDF controllers of the liquid slosh system. A considered amount of literature has been published on 
PSO tune PID in a number of applications such as [22, 23]. Those studies have proved that PSO is a good 
choice to tune the PID controller in their application. A nonlinear liquid slosh model in [24] which consists of 
a small motor-driven liquid tank performing a rectilinear motion is considered. The PSO tuning method is 
compared by heuristic tuning method in order to show the effectiveness of the proposed tuning approach. 
Based on the simulation results, the suggested tuning method is capable to reduce the liquid slosh level in the 
same time produces fast input tracking of the tank without precisely model the chaotic motion of the fluid. 

The paper is organized as follows: Section 2 described the research method which consist of 
modelling of liquid slosh, PIDF controller and PSO algorithm, Section 3 presents the simulation results and 
discussion and finally, Section 4 concludes this paper. 


2. RESEARCH METHOD 
2.1. Modelling of Liquid Slosh 

Figure 1 shows the modelling of liquid slosh that perform a movement in horizontal axis as a basis 
for simulation works in developing the tuning method for PIDF controller. The spring-mass-damper and a 
pendulum basically used to model the liquid slosh as a physical interpretation of fluid surface then created a 
slosh. The liquid slosh is modelled by a simple pendulum system performing the horizontal motion with a 
length, l and slosh mass, m. The slosh angle, @ is indicated by the pendulum angle as shown in Figure 2. 


By referring to Figure 2, 
Y = lsin + y (1) 
Z = l — lcos0 (2) 


The kinetic energy is defined as 
eee i GOR. A 2 
T =>M(y)* +>m(Y) +=m/(Z) (3) 
The potential energy is formulated as 
V = —mglcosé (4) 
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The Lagrangian of the system is the difference of kinetic energy and potential energy, 
L=T-V (5) 


Then, the Euler-Lagrange equations in y and 0, which produce dynamic equations of the system, 
is given by 


My + mlcos66 — ml6?sin@ = u (6) 
mlcos0y + ml?6 + dð + mglsiné = 0 (7) 


where d,g,0,Z,Y,y,u,l,mand M are damping coefficient, gravity, pendulum angle (slosh angle), 
displacement of m in the vertical direction, displacement of m in the horizontal direction, displacement of 
rigid tank, force applied for translational motion, hypotenuse length of the slosh (length of pendulum), 
mass of pendulum (slosh mass) and mass of the tank and liquid, respectively. 

Table 1 shows the value of the liquid slosh system parameters used in this study. The quick-stop 
experiment as mentioned in [25] is used to find the system parameters by considering the liquid 
characteristics, tank dimension and the ratio of liquid fill. 





„Rigid 
Tank 
Cart 
ue > 
b—>y 
Figure 1. Liquid slosh model Figure 2. Liquid slosh modeled by simple pendulum 


Table 1. Liquid Slosh Parameters [25] 


Parameter Value Unit 
M 6.0 kg 
m 1.32 kg 
l 0.052126 m 
g 9.81 ms~? 
d 3.0490 x 107 kgm?/s 


2.2. PIDF Controller 
Figure 3 shows the block diagram of liquid slosh problem with two PIDF controllers, one for liquid 
slosh angle and another one for cart position control where O(t), y(t), u(t), e(t) and r(t) are slosh angle, 
cart position, control input, error and reference, respectively. G, is the motor-driven liquid tank system. 
The transfer function of PIDF controller, K;(s)(i = 1,2) is defined as follow 
Ni 


K,(s) = P; +1,=+Di— (8) 


where P; is the proportional gain, I; is the integral time, D; is the derivative time and N; is the filter 
coefficient. 

The tuning method for PIDF parameters of liquid slosh control problem is done by minimizing the 
Sum Square Error (SSE) and Sum Absolute Error (SAE) based on the liquid slosh angle and cart 
position response. 


Indonesian J Elec Eng & Comp Sci, Vol. 18, No. 2, May 2020 : 750 - 758 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 g 753 





Figure 3. Block diagram of liquid slosh problem with PIDF controller 


2.3. Particle Swarm Optimization (PSO) 

PSO is a Swarm Intelligence (SI) algorithm motivated by the swarm of birds and fish schooling. 
In PSO algorithm, particles as known as agents will travel to it current best location and global best randomly 
by its current velocity that will be updated randomly later. Hence, to increase the efficiency of the swarm, 
the global communication between particles and real-number randomness is used in the algorithm. 

Let v; and x; is the velocity vector and position vector for agent or particle i, respectively. 8 PIDF 
parameters to be adjusted, so the design parameter is stated in (9) 


x; = [P, I, Dı Ny Pz Iz Dz N3] (9) 


Firstly, the particles velocity, vjand swarm particles, x; are set randomly using (10) and (11), 
respectively, where up is an upper boundary and l, is a lower boundary set for the agents, while 7, and r, is 
random real-number between [0 1]. 


Xi = lp Gy 1g) Kx (10) 
=F (11) 


For this simulation, an enhanced PSO algorithm using a priority-based fitness technique is 
recommended for PIDF tuning [26]. Sum squared error (SSE) of slosh angle, 0 is set as uppermost primacy 
followed by sum absolute error (SAE) of 0, sum square error (SSE) of cart position, y, sum absolute error 
(SAE) of y, overshoot (OS) of y and settling time, Ts of y. The main concern is to get the smallest number of 
these fitness functions which are as shown in (12) to (17). 


f, (x) = Fitness, = SSEg = } eọ° (t)dt (12) 
fa (x) = Fitness, = SAE, = dileg|(t)dt (13) 
fa (x) = Fitness; = SSE, = Yiey*(t)dt (14) 
fax) = Fitness, = SAE, = dle, |(Hat (15) 
fs (x) = Fitness; = OS, = Ymax — Yfinal (16) 
fo (x) = Fitness, = Settling Time = Ts, (17) 


The personal best agent is updated based on (18) as 


xi, fi (xt) < fA (XE) AND fa (xi) < fa (X£) AND... 
wo fg(xi) < FXE) AND falt) < GXH AND 
e fa (xi) < XD AND folt) < f6XD 


x;, otherwise 


(18) 


* == 
Xt+1 7 


The global best is updated by using in (19) 
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xi, fi (xi) < FCD AND fa (xi) < fa (gi) AND..... 
fai) < BOD AND falxt) < fagi) AND... (19) 
ae fs(xt) < fs Cgi) AND f(x!) < f(O 


gč, otherwise 


* E 
Jt+1 = 


Then, the updated velocity for the particle is obtained by (20) 


vitt = wX vita xr x (g*— xt) +B xn x (x* — xt) (20) 
where, an inertia function, w; as (16) is used to update the velocity, vj** in every new iteration (t + 1) per 
(22) [27], 7,3 and 7, is random real-number between [0 1], the parameters 6 and a are the personal coefficient 
and social coefficient, respectively. Typically, these coefficients are constant and set to 2. 





w; = 0.4 + 0.5 x 2 (21) 
The new position then is updated by 
ae ee) (22) 


The simulation stopped after the maximum iteration is reached and the best PIDF parameters 
obtained by PSO algorithm is presented. 


3. RESULTS AND ANALYSIS 

The performance evaluations of the proposed tuning method are based on the ability of the cart to 
follow the input in horizontal motion and liquid slosh level reduction in time domain. The step input is set to 
0.5 meter with zero initial conditions as stated in (23) 


0, O<t<0O.5, 

ee Gece (aa 

The simulation works are done with MATLAB 2016, Microsoft Window 10, 8GB RAM and Intel 
Core 17-6700 Processor (3.41GHz). The liquid slosh model as stated in (6) and (7) are designed via Simulink. 

The simulation is run for 20s and the sampling time is fixed at 0.01s. The maximum iteration for this 
work is set at only 12 iterations. For inertia function, w, the initial value is set to 0.9 and decreased linearly to 
0.4 by using (21). Only 2 particles or agents are used in this work that contribute to 24 number of 
evaluations. The upper and lower boundary is set initially as follow: 


lower boundary = l, = [0 0.001 0 0 100 50 20 100] 
upper boundary = u, = [5 0.005 10 5 120 60 30 120] 


The PSO tuning method is compared by manual heuristic tuning method in order to show the 
effectiveness of the proposed tuning approach and the results for optimal PIDF paratmeters for both methods 
is shown in Table 2. 


Table 2. PIDF Parameters Obtained by PSO and Heuristic Tuning 
PIDF gain PSO Tuning Heuristic Tuning 


P, 3.2769 20.0 
1, 0.0015 0.1 

D; 6.2659 10.0 
N; 1.3856 10.0 
P, 105.4218 110.0 
I 53.3039 60.0 
D, 20.3509 30.0 
N, 114.0931 120.0 
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Figure 4 shows the simulation result on cart position response for both PIDF controller tuning 
method. It can be seen from the response in Figure 4 that, for PSO tuning, the cart is able to follow the 
prescribed input trajectory presicely without any steady state error and small overshoot about 3%. 
The settling time recorded is 7.5 s. On the other hand, for the heuristic tuning, the settling time is much 
longer about 10 s and the response produces huge overshoot about 20%. In term of rise time, heuristic tuning 
produces faster rise time which is 1.5 s compared to PSO tuning which is 3s. What is interesting about the 
result in Figure 4 is PSO tuning produces much better performance compared to heuristic tuning in terms of 
time response specification. 

Figure 5 presents the slosh angle response for both tuning method. As can be seen from Figure 5, 
by heuristic tuning, the maximum residual of liquid slosh is +0.11 radian and settled within 9 s. Interestingly, 
by PSO tuning, there is significant improvement for liquid slosh reduction where the maximum residual is 
decreased to +0.025 radian and settled within 6 s. This result is supported with the slosh rate in Figure 6 
where by using PSO tuning, the maximum slosh rate is +0.25 radian/sec while by tuning heuristically, 
the maximum slosh rate is +2.3 radian/sec. The result of the necessary control efforts is shown in Figure 7. 
From the result we can see that the control signal overshoots when the desired stimulus is applied to the 
system. The maximum control input of heuristic tuning is bigger compared to PSO tuning. For heuristic 
tuning, the maximum control input is 60 Newton while for PSO tuning, the maximum control input is 5 
Newton. In summary, these results show that the PSO tuning for PIDF parameters is surpassed the heuristic 
tuning in suppressing the liquid slosh in the same time, achieved the desired cart position. Table 3 
summarizes the results obtained from the simulation analysis. 

















Figure 4. The response of cart position Figure 5. The response of slosh angle 





Shbsh (ato tadam eer 








Figure 6. The response of slosh rate Figure 7. The response of control input 


Table 3. The results of Liquid Slosh System 


Method of Maximum Rise Time, Maximum Steady State Control input Percentage Settling 
tuning Slosh angle, Tr (s) Slosh rate error Overshoot Overshoot, Time, Ts (s) 
0 (radian) (rad/sec) (N) POS (%) 
PSO 0.025 3.0 0.25 0.0 5 3 iPS, 
Heuristic 0.11 1.5 2.3 0.0 60 20 10.0 
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4. CONCLUSION 

The aim of the present research was to develop the tuning method for PIDF controller for liquid 
slosh system by implementing the PSO algorithm. This study has shown that there is a significant 
improvement of PSO tuned PIDF controller if compared to tuning heuristically in terms of cart capability to 
follow the desired input and reduction of liquid slosh in the container. The results of this study indicate that 
the proposed tuning method for PIDF controller is a good choice to be apply in aforesaid sectors in order to 
solve the liquid slosh problem. The tuning method used for this study may be applied to other nonlinear 
applications elsewhere in the world. The main weakness of this study was just a few number of particles and 
iterations is applied for tuning process. This is because the highly nonlinear model of liquid slosh that 
contribute to error during simulation works. More research is required to determine the efficacy of PSO 
tuning especially in the high number of particles and iterations. Further works on experimentation of liquid 
slosh tank driven by motor by using model-free control approach is strongly recommended. 
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